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The invisible signals of right-handed sneutrino decays originating from a Z' are analyzed at the 
Large Hadron Collider. The possibility of accessing these events helps disentangling the B — L 
extension of Minimal Supersymmetric Standard Model from more popular scenarios of Supersym¬ 
metry. We assess the scope of the CERN machine in establishing the aforementioned signatures 
when accompanied by mono-jet, single-photon or Z-radiation probes through sophisticated signal- 
to-background simulations carried out in presence of parton shower, hadronisation as well as detector 
effects. We find substantial sensitivity to all such signals for standard luminosities at Run 2. 


The B — L Supersymmetric Standard Model (BLSSM) 
represents an appealing non-minimal realization of Su¬ 
persymmetry (SUSY), as it is more compatible with cur¬ 
rent Large Hadron Collider (LHC) data than the Minimal 
Supersymmetric Standard Model (MSSM) and further 
uniquely accounts for the well established existence of 
non-zero neutrino masses. In this scenario, which is based 
on the gauge group SU(3)c xSU(2)l xU(1)y xU(1)b-l, 
(heavy) right-handed neutrino Superfields are introduced 
in order to implement a Type I seesaw [1], Also, it has 
been shown that the scale of the B — L symmetry break¬ 
ing is related to the soft SUSY scale [2]. Furthermore, 
the BLSSM also alleviates the little hierarchy problem 
of the MSSM, as both the additional singlet Higgs state 
and right-handed (s)neutrinos release additional param¬ 
eter space from the LEP, Tevatron and LHC constraints. 

The Z' and (s)neutrino sectors are ideal hallmark man¬ 
ifestations of the BLSSM as candidate underlying SUSY 
model. An intriguing signal would be the one involving 
totally invisible decays of a Z' into (s)neutrinos, thereby 
being potentially accessible in mono-jet, single-photon 
and Z-ISR (Initial State Radiation) analyses. Contrary 
to SUSY models which do not have a Z' in their spectra, 
in the BLSSM one can afford resonant Z' (with a TeV 
scale mass) production and decay into heavy (s)neutrinos 
which can in turn decay, again on-shell, into an invisi¬ 
ble final state. Therefore, one would expect the typical 
BLSSM distributions of the visible probe (whether it be 
mono-jet, single-photon or Z-ISR) to be substantially dif¬ 
ferent from other SUSY scenarios where the mediator is 
much lighter. In this letter, we aim at studying this phe¬ 
nomenology and assess the scope of the LHC in testing 
it. It is worth mentioning that visible decays of the Z' 
into (s)neutrinos at the LHC have been studied in [3]. 

The particle content of the BLSSM includes, in addi¬ 
tion to the MSSM fields, three chiral right-handed Su¬ 
perfields (Ni), a vector Superfield associated to U(1)b-l 
(Z') and two chiral Higgs singlet Superfields (xi, X 2 )- 
After B — L symmetry breaking has taken place, the 
U(1)b-l gauge boson acquires a mass [1], M§, = 


g 2 B _ L v ’ 2 , which is function of the B — L gauge coupling 
and the Higgs singlet vacuum expectation value. The 
most stringent constraint on the U(1)b-l gauge boson 
parameters is obtained from LEP2 results, which imply 
Mz ' > 6 TeV [41. However, one should note that this 
result is based on the assumption that the Z' dominantly 
decays to the Standard Model (SM) quarks and leptons 
and is derived from the limit on the low energy four- 
fermion contact interactions induced by the Z' exchange 
diagram. As such, this type of bound cannot per se be 
relaxed, as one cannot account for, e.g., invisible decays 
of the Z' (into right-handed (s)neutrinos, in the BLSSM). 
In fact, if one departs from the effective Lagrangian ap¬ 
proach, other dynamics may come into play, e.g., non- 
negligible interference effects between the Z' and the SM 
neutral gauge bosons (again, as typical in the BLSSM). 
All this being said, since a re-computation of the LEP 
limits in the BLSSM is beyond the scope of this work, 
we adopt the 6 TeV limit as conservative in our scenario 
and nonetheless use it to constrain the gs-L coupling [5], 
which should then be <?b-l < Mz >/6 TeV. 

The Z' parameters should also be consistent with the 
exclusion limits imposed by the LHC Run 1. The lat¬ 
ter are obtained from the Drell-Yan neutral channel as 
a function of the dilepton Branching Ratio (BR) of the 
Z ', which in the BLSSM depends on a variety of pa¬ 
rameters. However, rather than scanning on the latter 
(which will be the scope of a separate publication [6]) 
we have selected here two representative configurations 
of the BLSSM, one with a narrow and the other with 
a wide resonance, namely, Tz>/Mz> — 2% and ~ 40%, 
respectively. In case of a narrow resonance, one finds 
that Mz> > 2.5 TeV whereas, for the wider resonance, 
this bound can be lessened to 0(1 TeV), both for gs-L 
couplings that we have verified to be not excluded by 
the data collected at the LHC with a center of mass en¬ 
ergy of 8 TeV and an integrated luminosity of 20 fb _1 . 
In order to perform such LHC tests, we have first re¬ 
produced experimental searches appropriate for both the 
scenarios of narrow and wide resonance in the dilepton 
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final states (with both electrons and muons). In both 
cases we have included Next-to-Next-to-Leading Order 
(NNLO) QCD corrections as well as the quoted accep¬ 
tances and efficiencies in [7] for the electron and muon 
channels. We have finally combined the signal and back¬ 
ground rates obtained from the two leptonic signatures 
and verified that our benchmarks have a significance be¬ 
low 2. In the narrow case we integrated the signal in a 
small mass region around the peak, chosen in order to 
maximize the respective significance, the latter obtained 
using Poisson statistics. In the wide case, since the peak 
structure in not visible, we have adopted a “counting ex¬ 
periment” approach. We have integrated the new physics 
signal imposing a low cut in the invariant mass distribu¬ 
tion and have verified that the excess of Z' events (above 
the SM ones) is not sufficient to claim a discovery, adopt¬ 
ing Gaussian statistics. Similarly, we did with regard to 
the bounds in [ 8 ]. In the remainder we will focus on 
a narrow Z', for sake of illustration. Further, for our 
chosen benchmark, all sparticle states, especially the col¬ 
ored ones, are heavy (but the sneutrinos and some of the 
charginos/neutralinos entering the Z' decay chains yield¬ 
ing invisible signals), so that no sizable intrinsic SUSY 
background exists. 

From the BLSSM Lagrangian, the relevant interactions 
for the right-handed sneutrinos are given by 

Ant = (Yv)ijhPR{Vk2Xk )^(A H )cr/FR a (1) 

+ ('A) ij {U MNs)ii Vi Pr (N^ xl)(T„ R )j a VR a 

+ (Y v )ij(M N )j cos P(T Ll ) pilpH + (T VR ) (xj v Ra . 

Here, the rotational matrices T r l and Fj, r are defined as 
Tl = (r LL ,r iH ) and T v = (r„ R ,r„ R ). From these in¬ 
teractions, it can be easily concluded that, if the lightest 
right-handed sneutrino is lighter than the lightest slepton 
and lightest cliargino, then it decays into light SM-like 
neutrinos and lightest neutralinos. This decay channel 
would be invisible, since both light neutrinos and light¬ 
est neutralinos would be escaping the detector. Hence, 
given the discussed SUSY construction peculiar to the 
BLSSM, this pattern can provide a robust signature of 
B — L sneutrinos through mono-jet, single-photon and/or 
Z-ISR, further considering that the high Z' mass will 
force the missing transverse energy (and hence the pr of 
the visible probe) to be harder than in Z (into neutrali¬ 
nos) mediated events typical of the MSSM. 

In our calculations we have used SARAH [9] and 
SPheno [10, 11] to build the BLSSM and calculate 
masses, couplings and Branching Rations (BRs). We also 
considered the following benchmark (given at the SUSY 
scale): M z > = 2.5 TeV, g B -L = 0.4, Mp Ri ~ M Vr2 ~ 
M Prs ~ 580 GeV, ~ M Prb ~ M ^ ~ 740 GeV, 
m-± ~ 4, 0.9 TeV, m^o ~ 440 GeV and slepton masses 
of order 700 GeV. Furthermore, the matrix-elements for 
the parton level signals and backgrounds were derived 
from MadGraph5 [12]. Then, for showering and hadro- 
nisation we have used PYTHIA [13] whereas we have 
performed the fast detector simulations with PGS4 [14]. 


As jet finding algorithm, we have used a cone (with 
calorimeter cluster finder) of size A R = 0.5. Then, we 
have manipulated the generated data with MadAnaly- 
sis5 [15]. Finally, we have adopted usual selection strate¬ 
gies, wherein cuts are enforced against the kinematics of 
the highest pr jet/photon or the Z. 

We start with the mono-jet case, which is generally 
dominated by sneutrino decays [25] (hereafter, j =jet), 

qq ->• Z'{-+ UrUr -> XiXi vv) + j■ (2) 

The SM backgrounds are the following: Z(—> vv) + j 
(irreducible) plus W(—> Iv) + j, W(—> tv) + j, tt and 
ZZ(—> 2v2v) + j (all reducible). We closely follow 
here the selection of [16, 17]. Further, in order to in¬ 
crease the Monte Carlo efficiency and thus obtain rea¬ 
sonable statistics, we have applied a parton level cut of 
PtUi) > 120 GeV for both signals and backgrounds (here 
Pr(ji) is the highest jet transverse momentum). Accord¬ 
ing to the estimation of the QCD background based on 
the full detector simulation of [18, 19], such a noise can 
be reduced to a negligible level by requiring a large Ar 
cut. Thanks to the heavy Z' mediation, we can afford to 
set here At > 500 GeV [20]. In view of this we can then 
also implement a $ T >100 GeV cut for both signals and 
backgrounds at generation level. The beneficial effect of 
the At selection is evident from the left plot in Fig. 1 . In 
contrast, a similar cut on pr(ji) is not as selective (and 
is anyway correlated), see the right plot in Fig. 1, yet it 
pays off to also enforce it. 

In Tab. I, we present the actual cut flow for signal and 
background events, given at 14 TeV with an integrated 
luminosity of 300 fb -1 . After the BLSSM specific cuts, 
i.e., Ar and pr(ji) > 500 GeV, all the backgrounds are 
reduced under the dominant sneutrino signal. Other key 
steps of the analysis are the lepton, r- and 6 -jet vetoes, 
which suppress the W + j and tt backgrounds by more 
than two orders of magnitude [ 21 ]. 

Now we turn to the single-photon signature, which oc¬ 
curs in our BLSSM point mainly via the following pro¬ 
cess: 

qq -> Z'{-+ VrV* r -> XiXi^) + 7- (3) 

We have generated single-photon events while requir¬ 
ing the following parton level (generation) cuts: $ T > 
50 GeV, prill) > 40 GeV, p T {ji) > 25 GeV (prim) 
being the leading photon transverse momentum). We 
also generate the background processes Z(—» vv) + 7 
and W(—> Ivi) + 7 , where l = e, p or r. The latter 
are reduced significantly by applying the cut flow shown 
in Tab. II, see [16, 22] for guidance. Again, the sneu¬ 
trino signal emerges dramatically over the backgrounds. 
Further, in Fig. 2, the spectra in $ T and prill) are 
shown. These plots well motivate our high transverse en¬ 
ergy/momentum cuts, however, we note that the smaller 
signal rates here force us to a softer, yet still very effec¬ 
tive, requirement, of a 150 GeV threshold. 

For the Z-ISR signature in the BLSSM, 

qq —> Z'{—¥ vrVr —t XiXi u v) + Z , (4) 
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FIG. 1: (Left panel) Number of events versus the missing transverse energy. (Right panel) Number of events versus the 
transverse momentum of the leading jet. Distributions are for the mono-jet case given after the jet selection only. The energy 
is 14 TeV whereas the integrated luminosity is 300 fb -1 . Here, M Z ' = 2.5 TeV and gB-L = 0.4. 



Backgrounds 

Signals 

Process 

Z{vv)j 


W(rv T )j 

tt 

ZZj 

Z’(2x2u)j 

Z'(uD)j 

Z'(XX)3 

Before cuts 

21573000 

19248000 

9390000 

179058000 

6621 

1334400 

278 

7.54 

+3 

3 

O 

(I) 

16823567 ± 1924 

15817945 ±1678 

7719914 ±1171 

151390826 ±4836 

5732 ± 28 

1219314 ±324 

255 ± 4.68 

6.895 ±0.77 

( 2 ) 

65275 ± 255 

135191 ± 366 

65423 =h 254 

298430 =h 545 

73 d= 8.5 

130636 ± 343 

27 ±4.95 

0.741 ±0.82 

(3) 

45530 ± 213 

32569 ± 180 

27102 ±164 

6836.8 ±82.7 

55.6 ±7.43 

118456 ± 328 

25 ±4.74 

0.672 ±0.78 

(4) 

14283 ± 119 

10566 ± 102 

8668.5 ±93.1 

2808 ± 53 

16.5 ±4.06 

35424 ± 185 

7.4 ±2.68 

0.201 ±0.44 

(S) 

10831 ± 104 

7395.3 ± 86 

6088.7 ±78 

881.7± 29.7 

12.2 d= 3.49 

23330 ± 151 

4.9 ±2.18 

O.132±0.36 

(6) 

8992.5 ±94.8 

6o07.4±77.5 

4699.9 ±68.5 

379.8 d= 19.5 

9.79 d= 3.13 

18806 ± 136 

3.9 d= 1.96 

0.107 ± 0.33 

(7) 

8969.8 ±94.7 

3343.1 ±57.8 

3929 ±62.7 

257.7 ±16.1 

9.78 ± 3.12 

18786 ±136 

3.9 ± 1.96 

0.107 ± 0.32 

(8) 

8969.8 ±94.7 

871.2 ±29.5 

3207.4 ± 56.6 

176.3 d= 13.3 

9.77 ±3.12 

18782 ±136 

3.9 d= 1.96 

0.107 ± 0.32 

(9) 

8458.9 ±92 

790.2 ±28.1 

1378.8 ±37.1 

81.39 ±9.02 

9.21 ±3.03 

17878 ± 132 

3.7 ± 1.92 

0.102 ± 0.32 

(lu) 

8152.3 d= 90.3 

769.9 ±27.7 

1334.4 zh 36.5 

54.26 ± 7.37 

8.8 ±2.96 

17357 ±130 

3.6 =h 1.89 

0.098 ±0.31 


TABLE I: The cut flow on background versus signal events for M Z i = 2.5 TeV and gB-L = 0.4 in the mono-jet channel 
at y/s = 14 TeV with Cdt = 300 fb -1 : (1) n(jets) > 1 with 1 77 (ji)| < 2; (2) pr(ji) > 500 GeV; (3) $ T > 500 GeV; (4) 
A(fr(j 2 ,$ T ) > 0.5; (5) veto on prCte) > 100 GeV, \r]{j 2 )\ < 2; ( 6 ) veto on prOh) > 30 GeV, \r](j 3 )\ < 4.5; (7) veto on e; ( 8 ) 
veto on g; (9) veto on r-jets; (10) veto on 5-jets. 


we generate events with the following parton level cuts: 
$ T > 80 GeV, p T (l) > 10 GeV and p T {j) > 20 GeV. 
The dominant irreducible background is ZZ —> l + l~uu 
(l = e, n) and the other noise in this category is WW —>■ 
l + vl~v. As we reconstruct the Z probe, we enforce a cut 
on an invariant mass window centered on the Z mass for 
two oppositely charged leptons, mu £ [76,106] GeV [23], 
the latter is strongly reduced. The reducible backgrounds 
may have jets produced: Z + jets, ZZ —» qql + l~ and 
ZW —» l + l~qq. After the customary large cut that 
the heavy Z' allows us to enforce (here, > 250 GeV), 
this cumulative noise yields no event for the luminosity 
adopted. In addition, there are other reducible back¬ 
grounds with jets that we have dealt with: tt —1 l + vbl~Vb 
which is reduced by rejecting events containing at least 
one jet with pt > 25 GeV; FF+jets which is reduced 
by the large $ T cut. The last leptonic background 
is ZW — > lvl + l~ , which is also eliminated by the cut 
$ T > 250 GeV. The cut flow (modeled on [24]) and indi¬ 
vidual responses of signals and backgrounds are shown in 
Tab. III. In Fig. 3, we show the various distributions in 


$ T , again, for the purpose of justifying our BLSSM spe¬ 
cific ]fc T selection. Again, the invisible signal dominated 
by the sneutrinos and accompanied by the reconstructed 
Z stands well above the backgrounds. 

In summary, we have proven the sensitivity that the 
LHC has in Run 2 with standard luminosity settings in 
probing invisible signals which emerge in the BLSSM 
from Z' decays in presence of an associated jet, pho¬ 
ton or Z-boson. For all such signatures, we were able, 
upon enforcing well established selection procedures for 
these topologies supplemented by BLSSM specific cuts, 
to establish signals with significances well above the cus¬ 
tomary 5cr discovery limit. Indeed, this has been possible 
thanks to the fact that the BLSSM mediator of such in¬ 
visible signals is a very heavy Z\ with mass in the TeV 
region, thereby transferring to its decay products large 
transverse momenta that can be generically exploited in 
all cases for background reduction. Furthermore, for all 
topologies considered, the dominant component of the 
signal is via sneutrinos (above neutrinos and neutrali- 
nos), so that assessing these invisible signatures in the 
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FIG. 2: (Left panel) Number of events versus the missing transverse energy. (Right panel) Number of events versus the 
transverse momentum of the leading photon. Distributions are for the single-photon case given after the jet selection only. The 
energy is 14 TeV whereas the integrated luminosity is 300 fb _1 . Here, M Z ' = 2.5 TeV and gB-L = 0.4. 



Backgrounds 

Signals 

Process 

Z(yv)^f 

W(lvi) 7 

Z'{ 2*2i/)7 

Z'(yV) 7 

Z'(XX) 7 

Before cuts 

332712 

204644 

37380 

0.861 

0.234 

Cut 

n(y) > 1 

316031 ±125 

192677 ±106 

34998 ± 47.2 

0.806 ±0.227 

0.219 ± 0.118 

Pt( 7 i) > 150 GeV 

18576 ±132 

12146 ±106 

12357.8 ±91 

0.282 ±0.435 

0.0765± 0.2268 

$ T > 150 GeV 

14681±118 

4287.3 ±64.8 

11202 ± 88.6 

0.255 ± 0.424 

0.0693 ±0.2208 

n(j) < 1 , \r)(j)\ < 4 -5 

6819.7 ±81.7 

2388.2 ±48.6 

7415 ±77.1 

0.168 ±0.368 

0.0457 ±0.1917 

veto on e 

6817.6 ±81.7 

1731.8 ±41.4 

7409.3 ±77.1 

0.168 ± 0.368 

0.0456 ±0.1916 

veto on fi 

6817.6 ±81.7 

1132.5 zb 33.5 

7407 ±77.1 

0.168 ± 0.368 

0.0456 ±0.1916 

veto on r-jets 

6479.8 ±79.7 

758 ± 27.5 

7069 ±75.7 

U.151 ± U.531 

0.0435 ± 0.1882 


TABLE II: The cut flow on background versus signal events for M Z ‘ = 2.5 TeV and gB-L = 0.4 in the single-photon channel 
at y/s = 14 TeV with Cdt = 300 fb _1 . 


heavy ]fc T regime would not only signal the presence of a 
dark matter induced channel within SUSY but also be a 
circumstantial evidence of a theoretically well motivated 
non-minimal version of it, the BLSSM. In a forthcoming 
publication [6], we shall show that the case made here 
for illustrative purposes using a benchmark with a 2.5 
TeV and rather narrow Z' can be extended to a large 
BLSSM parameter space volume (also covering lighter 
and/or wider Z's). 
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FIG. 3: Number of events versus the missing transverse en¬ 
ergy. Distributions are for the Z-ISR case given after the jet 
selection only. The energy is 14 TeV whereas the integrated 
luminosity is 300 fb -1 . Here, M Z i = 2.5 TeV and gB-L = 0.4. 
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